t—J model of coupled Cu 2 Os ladders in Sru-zCa^Ci^C^i 
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Starting from the proper charge transfer model for CU2O5 coupled ladders in Sri4_ 2! Ca a; Cu2404i 
we derive the low energy Hamiltonian for this system. It occurs that the widely used ladder t—J 
model is not sufficient and has to be supplemented by the Coulomb repulsion term between holes in 
the neighboring ladders. Furthermore, we show how a simple mean-field solution of the derived t—J 
model may explain the onset of the charge density wave with the odd period in Sri4_ a; Ca a; Cu2404i. 
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I. INTRODUCTION 

It is widely assumed that the two-dimensional (2D) t- 
J modelă is the correct model to describe the low-energy 
physics of the CuC>2 planes.— Consequently, many au- 
thors believe that the high-tcmpcraturc supcrconductiv- 
ity in the cupratcs can be cxplaincd by this model and 
it is merely the computationally challenging character of 
the model which leads to the lack of the understanding 
of the superconducting ground state (see e.g. Refs. 0). 

Similarly, it has been suggested that the t-J model 
dcfined on the ladder (called ladder t-J model in what 
follows) is the right model to describe the low energy 
physics relevant for the CU2O5 coupled ladder plancs 
of Sr 1 4_ 2 .Ca a; Cu2404i (SCCO)££ This is a very attrac- 
tive theoretical idea because: (i) the ladder t-J model is 
much easier to solve than its 2D counterpart and it has a 
superconducting ground state for some specific range of 
parametersji (ii) a superconducting ground state (under 
pressure of 3 GPa) was foundâ in the ladder planes of 
SCCO for x = 13.6. This may suggest that indeed the t- 
J model contains the essential physics needed to explain 
the superconductivity, at least in the ladders. 

In this paper we would like to question the above point 
of view. As we show below, the ladder t-J model is too 
oversimplified and thus not sufficient to describe the low 
energy physics of the ladder planes in SCCO. Actually, 
this can alrcady be inferred by comparing the experimen- 
tal observations in SCCO with the theoretical predictions 
for the ladders^ (i) a charge density wave (CDW) ground 
state with period 3 and 5 was observed, but (ii) 110 CDW 
state with even period has been found^ whereas (iii) the 
ladder t-J model may have a CDW ground state only 
with an even period i 10 ' 11 Therefore, we investigate this 
problem here by a systematic derivation of the proper 
t-J model for the coupled ladder system, extended for 
topological reasons by the interladder repulsive term, and 
discuss a simple solution of this model. 

The paper is organized as follows. We introduce the 
charge transfer Hamiltonian in Sec. UH Ncxt, in Sec. 



IIIII wc derive the low energy t-J Hamiltonian which con- 
tains the kinetic energy and the superexchange, similar 
to the ladder t-J model, and the intraladder and inter- 
ladder repulsion terms. In Sec. IIVI we examine the role 
of the Coulomb intersite repulsion. Finally, we present 
a numerical solution of the model in Sec. [V] and draw 
conclusions in Sec. IVII The paper is supplemented by 
two appendices where some of the mathematical details 
of the model derivation arc discussed. 



II. THE CHARGE TRANSFER HAMILTONIAN 

As the starting point we choose the multiband charge 
transfer Hamiltonian introduced before for the CU2O5 
coupled ladder geometry in SCCOi^ The model in hole 
notation reads, 

n = n + n 1 +n 2 , (i) 

^0 = —tpd (^jagUiag ~~ ctcrV iota ~F < ^iaa X iaa 

± 4a + H.c.) + Aj^ (n lax + n i( 

ia 

(2) 

i ia 

Ux = u p Y, 

nia^nia^ +U p ^2 n ibt n ibi , (3) 
ia,^—x,y i 



^2 2?y) ^ ^ [j^iaxa^'iaycr "T" ^iaya^icxxa 

iota 

~t~ Up{\ 3//) ^ ^ (^icixa^iăya ^^iaya^iaxa^j • (4) 



The model (1) was adopted to the present ladder 
geometry^ from the charge transfer models introduced 
before for Cu02 planes^ and CUO3 chainsi 3 . in high 
temperature superconductors. The parameters are: the 
energy for oxygen 2p a (2p x or 2p y with creation opera- 
tors {xî } and {y\ aa }) orbital A (the so-called charge 
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FIG. 1: Schematic view of the coupled C112O5 (white/gray) 
ladders in SCCO: (a) orbitals in charge transfer model {]]); 
(b) intraladder (interladder) bonds in the effective extended 
t-J model, see Eq. ([5|, shown by solid (dashed) lines. 



transfer cncrgy mcasured with respect to the energy of 
3c! copper orbitals), the d-p hopping t p d between the 
nearest neighbor copper and oxygen sites, the on-site 
Coulomb rcpulsion U (U p ) on the copper (oxygen) sites, 
and X] = Jh/U p ~ 0.2 — a realistic value of Hund's ex- 
change on oxygen ionsi^ (for a complete set of realistic 
parameters see Sec. IIIIE|) . Besides, in principie the ac- 
tual electron energy at bridge orbital (rung position) of 
the ladder with creation operators {b\ a } is approximately 
10 % smaller than the one at other oxygen positionsJ^ 
However, it was showr>i£ that this difference does not 
have any important physical consequences and thereforc 
we will neglect it here. Note also that the phases of the 
{3d, 2p] orbitals were explicitly taken into account in the 
hopping clements cx t pc [ (for clarity the phases are shown 
only in Fig. [5]), the index a £ {R, L} denotes the right 
or left leg of the ladder (R = L and L = R) 7 ă = —a for 
(7 € {Î4}, and the upper (lower) sign stands for terms 
with a = L {a = R). 

The charge transfer model (1) includes seven or- 
bitals per CU2O5 ladder unit cell i (see Fig. [ŢJ: two 
Cu{2>d x 2_ y 2 = d) orbitals on the R/L leg, two 0(2p y = y) 
orbitals on the R/L leg, two 0(2p x = x) side orbitals on 
the R/L leg, and one 0(2p x = b) bridge orbital on the 
rung. Although it seems that the model is quasi one- 
dimensional (1D), the density operators fi{ axa and fii av(J 
stand for the oxygen hole densites in the ncighboring lad- 
ders and make it implicitly 2D as the interladder coupling 
couples the ladders, so the model extends over the entire 
CU2O5 plane. 



thus also called t—J—V model, 

H = Ht + Hj + Hy 1 + Hy 2 



(5) 



which are separatcly derived and discussed bclow. 

We begin with the superexchange term Hj which is 
the only important term in Eq. ([5]) at half-filling. In this 
case (i.e. with one hole per copper site) the charge trans- 
fer model can be easily reduced to the low energy 
Heisenberg model for spins 5 = 1/2 using the perturba- 
tion thcory to fourth order in t p 



.18 



+ [ SlR '■' Sii ~ I ™iR™iL ) • ( 6 ) 



Here, tilde in hi a implics that the hole double occupan- 
cies are excluded. The superexchange constant contains 
contributions due to charge excitations on copper sites 
and on the intermediate oxygen site for a Cu-O-Cu bond, 
and for finite U p rcads^ 



J 



'2t 2 * 



1 



U 2A + U V 



(7) 



One may wonder whether the geometry of coupled 
ladders could influence the above result. Indeed, there 
exists a 90° superexchange process between the holes 
on two ncighboring ladders. However, aceording to 
the Goodenough-Kanamori rulesi^ such a superexchange 
process is much weaker than the superexchange gener- 
ated by charge excitations along the 180° path in the 
single ladder and can be neglected. 



B. Zhang-Rice singlets for the ladder 

Whcn the ladder is away from half-filling, the pertur- 
bation theory gets complicated. Therefore, following the 
Zhang and Rice construction^ we first define a phase co- 
herent symmetric plaquette stat£& \Piaa) which is formed 
by the four oxygen orbitals surrounding the central cop- 
per site ia: 

\P iaa ) - \ (±x\ aa T bl - y\_ Xiaa + yj aa ) |0) , (8) 



III. THE EFFECTIVE t-J HAMILTONIAN 

A. The model and the superexchange 

In what follows we will derive the low-energy version 
of Hamiltonian ([1} which is valid in the so-called charge 
transfer regime U > A, i.e. for the typical values of 
model (HJ parameters: U ~ &tpd, A ~ 3t pc ;, and U p ~ 
3t p d, see Ref. [Ţ3 and Sec. IUI El bclow. The effective 
t-J Hamiltonian consists then of four terms, and may be 



TABLE I: Binding energy of the singlet and triplet state 
formed by the copper hole and a doped hole in one of the three 
various oxygen states: (i) symmetric plaquette state \Pi aa ), 
(ii) antisymmetric plaquette state |Ai ao -) orthogonal to Eq. 
([5]), and (iii) single oxygen orbital. Here t 1 = tp d / A ~ t pc i/3, 
*2 = tld/{U - A) ~ t pd /5, and t 3 = hU p /(A + U p ) ~ t pd /6 
for the typical charge transfer parameters.— 

\Piaa) \A iatT ) single oxygen 

singlet -8(ti + t 2 ) + 2t 3 -4ti + 2t & -2(ti + t 2 ) + 2t z 
triplet -4ti 
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where the upper (lower) sign stands for a — L (a = R). 
When a hole in this state forms a singlet state with the 
hole at the central copper site, it has a large negative 
binding energy of — 8(ti + £2) + 2t3, see caption of Table I 
for definition of t„ hoppings and for more details. Actu- 
ally, this binding energy is not only much larger than the 
individual cffective hopping terms (which is of the order 
of ti or t2, see Ref. 0) but it is also considerably larger 
than the binding energy of some othcr possiblc bound 
states, see Table I. Note that finite U p , not considered 
by Zhang and Rice^ results in finite £3 hopping but does 
not qualitatively change the large binding energy of a 
symmetric singlet state ([8]). 

The mere problem with the states defined by Eq. (|8]) 
is that they are not orthogonal. It can be checked that in 
the case of the laddcr gcomctry the following supcrposi- 
tion of the symmetric plaquette states forms a complete 
and orthogonal basis for a low-energy Hilbert subspace: 

•O) = ^Y, etkle ~ lk3 ( a k\P ja a) +Pk\Piă*)) , (9) 
jk 

where a k (p k ) = 2/V3 - 2 cos k ± 2/^5 - 2 cos k. Then 
the Zhang-Ricc (ZR) singlets for the ladder are: 

Hha) = 71 (ti^U - ^U4 t ) 10} ■ (10) 

Although the binding energy is slightly reduced after this 
orthogonalization, the change is not significant: If the 
energy splitting between the orthogonalized ZR singlets 
and triplets is defined as 16x 2 ^i (we consider a simplified 
case t\ = t2 and U p = 0), then \ ~ 1 — both in the 1D 
(X = 0.98) and in the 2D case (x = 0.96), see Ref. 0. 

Having shown that the ZR singlets in the single lad- 
der do not differ much from those which arise in the 2D 
cupratesj^i we can now safely apply all the arguments 
used in Ref. to derive the effective hopping of ZR sin- 
glets following from finite t p d- Thus, we obtain, 

H t = -tJ2{dlÂ ă , + (4X+l,aa + H.C.) } , (11) 

iota 

where once again di aa = di acr (l — ni aS ) is a fermion op- 
erator in the restricted space. While we do not show 
here the detailed expression for the effective hopping t 
of ZR singlets, note that it is considerably smaller than 
t p d (ca. 50%).— Note also that having two ZR singlets 
at the same site costs energy 4^2 + 2ti (if ta = 0, sec 
Ref. 0) and therefore we used the tilde operators above 
to exclude these local configurations of two ZR singlets. 

C. Intraladder repulsion Hv 1 

The Coulomb interaction on oxygen sites U p , neglected 
111 Ref. i plays a minor role in the stability of the ZR 
singlets (see e.g. finite £3 for finite U p in Table I), but this 
issue is more subtle<22 Actually, due to finite U p the two 




FIG. 2: (Color online) The artist's view of the intraladder 
repulsion between two neighboring ZR singlets. Large (small) 
arrows depict the hole spins for +1.0 (+0.25) charge — they 
stand for spins at copper sites and for the spins of doped holes 
delocalized over oxygen orbitals. Orbital phases are depicted 
by striped/white areas. 

neighboring nonorthogonal ZR singlets repel each other 
when two holes occupy a common oxygen site (see Figs. 
[2] and [3j) . Whereas the significance of the interladder re- 
pulsion is discussed in the next subsection, let us concen- 
trate first on the repulsion between the ZR singlets within 
a single ladder (see Fig. [2]), and calculate repulsion cx U p 
between two orthogonalized ZR singlets within the lad- 
der (tpsaj ^rplUiliphŞ, ^ja}- Let us note that the 'mixed 
terms' such as (4> S R,tprL\'Hi\''PhL,4'jL), which could apri- 
ori destroy the ZR singlets, fortunately turn out to bc 
much smaller than the respective binding energy. Using 
Eq. (fTUţ) , after a somewhat lengthy but straightforward 
calculation (for more details see Appendix [5} , one finds 
the following values for the intraladder interaction along 
the leg and the rung: 

(ipj a ,tpj + i ta \Hi\ipj + i, a ,ipj a ) = 0.027Up, (12) 
(lpja,1pjă\Hl\lpjă,1pja) = 0.026 t/p. (13) 

We have verified that the interaction between the sec- 
ond nearest neighbors is ca. 15 times smaller and can 
be safely neglected (the longer-range interaction is even 
smaller, cf. Appendix lA"ţ . 

Hencc, one finds that the interaction among the nearest 
neighbor ZR singlets is almost isotropic. Thus, we can 
write the effective Hamiltonian for the repulsion between 
ZR singlets (as shown in Fig. [2} 

Hvi = Vi( y^ J n ia h i+ i ta + ^2 n lR n lL ^j , (14) 

ia i 

where V\ ~ 0.027t/ p . Let us note that the ratio V\/U p 
is approximately 14% smaller than the naively estimated 
nonorthogonal value 1/32 = 0.03125. We also checked 
that dimcnsionality drives the following trend in the ra- 
tios Vi/U p : 0.023, 0.025, 0.027 (considered herc), and 
0.029, for a single rung, the 1D ladder, and the 

2D case, so V±/U p inercases with the increasing number 
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FIG. 3: (Color ordine) The artist's view of the interladder 
repulsion between two ZR singlets on two different (white 
and gray) ladders. Spins are depicted similarly as in Fig. [2] 

of neighborsi22 This follows because in lower dimensions 
the charge escapes easier from the orbitals b and y (pro- 
viding the dominating contribution) , while in the 2D case 
all the orbitals suffcr from the orthogonality problem. 

D. Interladder repulsion Hv 2 

Finally, we calculate the interladder repulsion be- 
tween the ZR singlets due to on-site repulsion U p 
in orbitals belonging to two neighboring ladders: 
(tpsa^r+^pl^l^h+^^i^joc) — a bar sign over ip de- 
notes the singlet formed on the neighboring laddcr. Be- 
sides, since the neighboring ladder is misaligned by a 
lattice constant 1/2 with respect to the one considered, 
we labei the ZR singlets on the neighboring ladder by 
j + 1/2 (for the copper-copper lattice constant equal to 
1). Next, using Eq. (jTUj) one finds after a somewhat te- 
dious but straightforward calculation (for more dctails 
see Appendix [Bj the following value for the interlad- 
der interaction between the closest sites belonging to the 
neighboring ladders (see Fig. [3]), 

(^*^j±i,*\U2$ j±h>gi rt ja ) =0.136 (1-5*7/2) tf p , (15) 

while othcr (neglected) longer-range repulsive terms are 
at least one order of magnitude smaller, cf. Appendix 
|B1 Thus, the repulsion between holes on the neighboring 
ladders reads: 

Hv 2 = V 2 ^2 («ia^+I,â + "«Ă-i,â) > ( 16 ) 

ia 

where hi a operator is related to the ZR singlets as before 
and V2 — 0.136(1 — 5rj/2)U p . We again neglected all 
spin-flip terms which arc small in comparison with the ZR 
binding energy and give zero when 'sandwiched' in the 
singlet states. Besidcs, the numerical prcfactor (cqual 
to 0.136) is here slightly enhanced with respect to the 
expected 1/8 = 0.125 value (unlike in the intraladder 
case). This is because a significant fraction of charge 
escapes from the b and y orbitals to the x orbitals due to 
the orthogonalization procedurc. 



TABLE II: Adopted values of the parameters of the charge 
transfer model (fi} from Ref. [Ţ3 (Jh from Ref. ITir i and the 
calculated values of the derived t-J-V model ([5]) in eV. 



char^ 


re transfer model 




t-J-V model 


t pd 


1.3 


t 


0.54 


A 


3.6 


,1 


0.24 


u P 


4.0 


W 


0.11 


U 


10.5 


Vi 


0.27 


Jh 


0.8 







E. Parameters of the effective model 

The calculated parameters of the effective t-J-V 
model ([!]) derived above are shown in Table HB This cal- 
culation is based on the cuprate charge transfer model 
parameters from Ref. [T?] (the value of Jh is takcn from 
Ref. IÎ4T ) which might be considered as the most widcly 
accepted choice of the cuprate parameters, cf. Refs. l2ll 
and I22I. Due to the same Cu-O distances in SCCO as in 
CuC>2 planes in e.g. La2Cu04, we can can adopt these 
parameters also to the present case. 

Let us note that on the one hand, it should bc empha- 
sized that the interladder coupling V2 is of the order of J 
for the realistic parameter a 14 i 17 and therefore cannot be 
neglected. On the other hand, the value of V\ is two and 
a half times smaller and therefore we suggest that, if nec- 
essary, this interaction could be skipped in the first-order 
calculations. 



IV. ROLE OF THE INTERSITE COULOMB 
REPULSION V pd 

One may wonder whether the intersite Coulomb repul- 
sion V p d in the charge transfer model could alone lead to 
a significant repulsion (i.e., of the order of the estimated 
value of V2) between the ZR singlets in the neighboring 
ladders. This term, which stands for the repulsion be- 
tween charges situated in the nearest neighbor copper 3d 
and oxygcn 2p orbitals, was neglected in Ref. |9| and in 
the above analysis, cf. Eq. (TT]). In fact, including this 
term may e.g. lead to a significant renormalization of the 
parameters of the 2D t-J modelă 

Indeed, one finds that the repulsion between a hole in 
a symmetric state \4>i a ) and a copper hole in state \dj a i) 
situated on the nearest neighbor sites in two neighboring 
ladders is of the order of (n ixcr + n lytJ )n djcrl Vpd ~ 0-5V p d- 
Although typically V p d is smaller than t p d, e.g. V p d ~ 1-2 
or even V p d < 1 eV^l this contribution might still 
be significant and, in principie, should not be entirely ne- 
glected. However, the key observation is that this term 
leads to roughly equally large energy cost if: (i) either 
the two ZR singlets are situated on the above mentioned 
sites i and j and repel each other due to V pd , or (ii) the 
two ZR singlets are situated far away from each other on 
two ladders and (due to V p d) merely feel the repulsion 
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FIG. 4: (Color online) The self-consistently calculated order parameter p for SCCO as obtained from trie effective t—J—V model 
([5]) [see Eq. (|17|l and text for more details] as a function of the interladder interaction Vi for: (a) filling n — 2/3 and period 
A = 3, (b) n = 3/4 and A = 4, (c) n = 4/5 and A = 5. Parameters: realistic values (see text) J = OAt and Vi = 0.2t (solid 
lines), and J = and Vi = (dashed lines). 



with the neighboring copper holc oii the neighboring lad- 
der. Hence, including the intersite repulsion V p d increases 
the total energy of the system but almost does not con- 
tributc to the energy diffcrence bctween the two above 
situations, measurcd by the value of the interladder re- 
pulsion Vi. It is only a small rcsidual repulsion due to 
the orthogonalization procedure, see Eq. ^ , which may 
change the value of Vz by a small fraction. This has been 
also confirmed by the results of Ref. HH, where V p d leads 
indeed to a very small repulsion between holes in the 2D 
t—J model (ca. 0.05t p d). 

V. NUMERICAL RESULTS 

Now we shall verify whether the derived t—J—V model 
Eq. ([5]) supports the CDW states observed in SCCO and 
to understand to what extent the interladder tcrm (|16j) 
influences the stability of the CDW state. As a thorough 
investigation is beyond the scope of this work and left 
for future studies, we solve the model ([5]) in the simplest 
possible way. Thus, we first introduce the Gutzwiller fac- 
tors, g t = (2 - 2n)/(2 - n) and g,j = 4/(2 - n) 2 , which 
rcnormalize the kinetic (gt) and interaction (gj) tcrms; 
for thcir justification see e.g. Refs. Here n denotes 
the average numbcr of d holes per site in the effective 
model ([5]), i.e., n = \^2 a (ni aa ) . Second, we use the mean- 
field approximation for the interaction terms. Next, we 
diagonalize the effective one-particle Hamiltonian intro- 
ducing the classical fields 

. j . _ J n — p for i/X 6Z 

(<W t W ~ | n+ _L_ p toi/x^z ' ^ U) 

wherc p < n is the CDW order parameter and A is the 
CDW period. Furthermore, ao .<îi_A l£M r) are dehned 

asinEq. jŢTJl but with i/X replaced by (i+l)/X [(i+2)/X 
for A = 5] — this assumption minimizes the classical 
energy cost of the interladder repulsion V2 . 

In Fig. 2] we show the CDW order parameter p cal- 
culated self-consistently as a function of the interlad- 



der interaction V2 for the three experimentally inter- 
esting doping levels>^£ n = 2/3 which corresponds to 
n/j = 2 — n = 4/3 holes per copper site in charge trans- 
fer model (fîj, n = 3/4 corresponding to rih = 5/4, and 
n = 4/5 corresponding to rih = 6/5. The results demon- 
strate that the interladder interaction plays indeed a cru- 
cial role in the stability of the CDW while the intraladder 
one is rather unimportant. 

We have found that the CDW state with period A = 3 
(A = 5) is stable for n = 2/3 (n = 4/5) for the rather 
realistic values of the parameters J = OAt, V\ = 0.2i and 
Vi = 0.5i — 0.9i. Please note, that: (i) we adopted here 
a somewhat smaller value of J = OAt which is closer 
to a typical value for cuprates^ and (ii) values of V2 
exceeding 0.5t can be obtained usingfor example the 
set of parameters suggested in Ref. [24|. This finding 
explains well the experimental results of Ref. 0. Besides, 
the CDW ordered state is also stable for period A = 4 
which was not observed^ We expect that the stability of 
the CDW state with this period is a shortcoming of the 
above simplified solution which does not capture well the 
frustration between two possible CDW patterns in the 
neighboring laddcrs which oceurs for period A = 4. - 



VI. CONCLUSIONS 

In summary, we derived the t—J model which describes 
the low energy physics of CU2O5 coupled ladders. Apart 
from the 'standard' superexchange oc J and kinetic en- 
ergy terms oc t, the model contains also the repulsion 
between the nearest neighbor holes in a laddcr oc V\ and 
in the two neighboring ladders oc Va , and hence is also re- 
ferred to as a t-J-V model [see Eq. ([5])]. We showed that 
the latter V2 term is roughly two and a half times larger 
than V\ and (contrary to Vi) cannot be skipped — in fact 
it is crucial to explain the onset of the odd period CDW 
state in SCCO. We emphasize that this particular extra 
term is restricted to the copper oxides in which oxygen 
is coordinated by three copper ions in the same plane. 
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Thcrcfore, it is not present in the CuC>2 planes^ or in 
Cu-O chaini^ in copper oxides, but (apart from the dis- 
cusscd SCCO case) could become relevant for the couplcd 
chains of SrCu02 (provided they are hole-doped). Fur- 
thermore, it is both a many-body term and of the order 
of J, contrary to various corrections to a 1D or 2D t—J 
modcl4a22 Besides, we also verified that the Coulomb in- 
tersite interacton V p d alone (not includcd in the presented 
derivation) cannot lead to a significant intcrladdcr repul- 
sion V 2 . 

The simple mean-field solutions of the t—J—V model 
derived here provides evidence in favor of the experi- 
mental observations of the onset of the odd CDW state 
in the ladder planes of the SCCO. This is further sup- 
ported by the recent density matrix renormalization 
group calculations^i where a ladder t—J model with the 
interladder coupling V2 (denoted as V± in Ref. H3) was 
studied: also there the CDW with odd period is stabi- 
lized due to the presence of the interladder interaction 
V 2 . 

Finally, let us note that the other hole-dopcd ladder 
compound Lai_ 2 ,Sr 2 ,Cu02.5 is also characterized by a 
large (but differcnt than the one discussed here) interlad- 
der couplingj^ Therefore, we argue that it is currcntly a 
challenge for the condensed matter community to search 
for a hole-doped ladder compound which could indeed be 
modelled by the ladder t—J Hamiltonian [given by Eqs. 
(|6|) and i.e., without additional intersite repulsion 

tcrmsl . 
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Appendix A: Derivation of the intraladder repulsion 
term oc Vi 

Here we show how to calculate the repulsion between 
orthogonalized ZR singlcts within the ladder due to the 
on-site interaction U p in Hi- Thus, one necds to deter- 
mine the following matrix elements: 

('(p sa ,tpra\ / Hl\tpha,1pj a ), {ipsa, 1prâ\Hl\tphâ, i>ja)- 

(Al) 

Let us note that the mixed terms such as for example 
(ip S R,tprL\'Hi\iphL,4'jL) vanish in the ZR singlet basis - 
they could a priori lead to the destruction of the ZR 
singlets, but fortunately they are much smaller than the 
respective binding energy. 

Intraladder repulsion along the leg. — First, we calcu- 
late the matrix elements of Hi between the orthogonal 



plaquctte states Eq. ([9]) along the leg: 

(4>saer:4>raă\Hl\4>ha9,4>jacr) 

I ■ \ \^ ik(h-r) iq(j-s) if(r-s) 
16 U P N 3l^ e e e 

kqf 

x I — (a k a q + p k f3 q - a k (3 q - Pk^qJ 



x (a q ^ s a k+s + /3 q ^f/3 k+ f - a q - f /3 k+f - Ş q ^fa k +fj 

■ k . q . q-f. k + f 1 \ 

sin — sm - sin — - — sm — h — 

2 2 2 2 16/ 

a k a q a q _fa k+ f + [3 k f3 q f3 q _ff3 k+ f^ |, (A2) 



and 

{4>sam 4>ras\Hl\4>haa, 4>jas) 
= -(<t>sa<r>(t>ra,<j\'rll\<t>ha<x,<t>iao), (A3) 

while the same spin elements are zero, 

{<t>saa,<f>raa\'rll\<f>hacT, ( f>jctcT) =0. (A4) 

One can evaluate numcrically the above expressions. It 
oceurs that the largest positive element is the nearest 
neighbor interaction 

(4>jacr,4>3 + l,as\'Hl\4'j+l,a9,4 l jacr) = 0.0544C/ p , (A5) 

while following Eq. (|A3|) the absolute value of the largest 
negative element, which corresponds to spin-flip nearest 
neighbor interaction, is the same. Furthermore, the sec- 
ond largest clement is the next nearest neighbor interac- 
tion and is over 20 times smaller, which means that it 
can be safely neglccted. 

Second, we calculate the matrix elements of Hi be- 
tween the nearest neighbor ZR singlcts, defined by Eq. 
(fTTiţ) . This introduces a factor 1/2 to the above esti- 
mations of the repulsion between orthogonal plaquette 
states: It is because there is a 50% probability to have 
opposite spins on a particular shared oxygen site oceu- 
pied by two holes from two different ZR singlcts. Note 
that the spin-flip-plaqucttc terms do not give any con- 
tribution to the repulsion between ZR singlcts, although 
they could in principie destabilize the ZR states them- 
selves. Fortunately, this is not possible since the binding 
energy of the ZR singlets is much larger. Thus altogether, 
we obtain for the repulsion along the same leg 



(lpj a ,1pj +lta \Hl\lpj+l, a ,i>j 



0.0272 U„. 



(A6) 



Intraladder repulsion along the rung. — Following a 
similar scheme, one can calculate the repulsion between 
ZR singlcts on different legs. One obtains the following 
matrix elements of Hi between the orthogonal plaquctte 
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states Eq. ((9]) on differcnt lcgs 



and 



aă\7~L\ | 



J_rr \ \ "* ik(h-r) iq(j-s) if(r-s) 



fe<3/ 



+ 



{ 16 i OLk ^ q + ^ k0lq ~ ak<Xq ~ ^ k ^ q ) 

(atq-fPk+f + P q - f a k+f - a q - f a k+f - (3 q -fP k+f ) 

( ■ k . q . q- f . k + f 1 \ 

sm — sm — sm sm 

V 2 2 2 2 16/ 

(a k q ai q -fPk+f + Ş k a q Ş q -fa k +/j } , (A7) 



and 



(4*sacr 



uă\7~Ll\4>hua, 4>jaâ} 
,<t>rââ\Ul\(t>hăv,<t>jot<r), (A8) 



and 



sacr ţ (firaa |^1 1 



0. 



(A9) 



Evaluating numerically the above expressions one obtains 
that the largest element is the nearest neighbor repulsion 
— this time between the orthogonal plaquette states on 
the same rung: 

(</>ja<T,<f>j,ăâ\'Hl\<t>j,ăă,<t>jauT) = 0.0529 U p , (A10) 

whilc the second largest element (the next nearest neigh- 
bor interaction) is over 15 times smaller and can be ne- 
glccted. 

Finally, following the same steps as those leading from 
(|A5j) to ()A6j) , we obtain the repulsion between the nearest 
neighbor ZR singlcts [dcfincd by Eq. (|T0|) ] along the same 
rung which is twice reduced: 

(^a,^ ă |Hl|^, S ,V^a> = 0.0265t/ p . (AII) 



Appendix B: Derivation of the interladder repulsion 
terni oc V2 

Here the task is to calculate the repulsion between two 
ZR singlcts centered at the ncighboring coppcr positions 
of two ladders (and thus sharing the same oxygen sites 
but not the p orbitals, see Fig. [3]) duc to the on-site 
repulsion on oxygen sites. However, again we will calcu- 
late the repulsion between arbitrarily located ZR singlets 
and only then we will show which elements are negligi- 
ble. Note that the plaquette states on two ladders are 
orthogonal to cach othcr although they still havc to bc 
orthogonalized for the same ladder (as in Appendix O . 
Explicitly one needs to calculate the following matrix el- 
ements: 



(Bl) 



(B2) 



Interladder repulsion between plaquettes with the same 
spin. — We calculate the matrix elements of H.2 between 
the orthogonal plaquette states Eq. (j9j) with the same 
spin but situated on different legs: 

(<?Wr, s +!,a<rl^2|<?Wi,S<r! <?W) = 

-(1 - 3t])U p -^ u k a q a q _ f a k+f 



16 



kqf 



| — sin q sin(g — /) + — sin k sin(fc + /) | 



x gifc('î-' , ) e i<î(j-s) e i/(r-s-i) 



while for the same legs we obtain 

(<f>raa, <Ps+^,aa\^\$ h+ ^ QrT , <j>jcta) = 

— (1 - 3r?)C/ p -^3 UkP q a q ~fPk+ f 



(B3) 



kqf 



| i sin q sin(g - /) + i sin k sin(fc + /) } 



x e ik{h-r) e iq(js) e if(r-s-±) _ 



(B4) 



As it might have been expected, it oceurs that the biggest 
term is the repulsion between orthogonal plaquette states 
with the same spin situated on the closest possible sites 
in the neighboring ladders (see Fig. [3]): 

(<?W, 0j±§,cU^2|0j±â,â<7) <?W) = 0.1355 (1 - 3ri)U P7 

(B5) 

and all other terms are of the order of 10~ 3 (1 — 3rf)U p 
and can be neglected. 

Interladder repulsion between plaquettes with opposite 
spin. — A very similar calculation as above, but per- 
formed for the orthogonal plaquette states Eq. ([5]) with 
opposite spins leads to the repulsion between orthogonal 
plaquette states with opposite spins and situated on the 
closest possible sites in the ncighboring ladders: 

(<t>jaa , 4>j± I ,a<7 1^2 1 4>j± 1 ! ^jaa) = 0. 1355 (1 - 2l]) U p , 

(B6) 

while again all other longer-range repulsive terms can be 
neglected. 

Finally, combining Eqs. (|B5|) - (|B6|) with the definition 
of the ZR singlet (fTCTţ) we obtain the value of the repulsion 
between the two ZR singlets (cf. similar discussion in Ap- 
pendix rAJ on the closest possible sites in the neighboring 
ladders to be 



(^,$ j± iJH 2 $ j±iiă ^ ja ) =0.1355 (l-5r//2)t7 p . 



(B7) 
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